function Y = naturalSelection(X)
%natural selection
    Xrs = sortrows(X, -3);
    m = ceil(size(X, 1) / 2);
    Y = Xrs(1: m, :);
end